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CORRECTED-INTERCEPT  CONTROL  OF  TORPEDO  MK48 
WITH  CENTROID  TRACKING* 


Introduction 

Thfe  acquisition  models  that  are  described  and  developed  in 
reference  (a)  are  based  on  the  assumption  that  tracking  either 
is  prevented  by  the  torpedo  noise  (complete  masking)  or  can 
proceed  without  degradation,  depending  on  the  positions  of  the 
torpedo  and  target  submarine  relative  to  the  tracking  submarine. 

In  the  present  report  we  modify  the  models  to  include  centroid 
tracking,  in  which  the  sonar  operator  tracks  a  composite  signal 
at  a  bearing  that  is  a  weighted  average  of  the  bearings  of  the 
torpedo  and  the  target  submarine.  The  weight  depends  on  the 
relative  intensities  of  the  two  signals  at  the  sonar  receiver. 

As  the  distance  of  the  torpedo  from  the  tracking  submarine 
increases  and  no  control  is  exercised,  the  composite  bearing  line 
shifts  gradually  from  a  position  near  the  torpedo  bearing  line 
to  a  position  near  the  target  bearing  line,  under  most  tracking 
conditions.  The  motion  of  the  composite  bearing  line  when  control 
is  exercised  depends  on  the  control  mode  and  the  procedure  for 
computing  the  changes  in  the  torpedo  course  angle.  In  this 
report  we  examine  the  motion  of  the  composite  bearing  line  and 
the  effects  of  this  motion  when  corrected-intercept  control  is 
exercised  and  the  corrections  to  the  torpedo  course  angle  are 
computed  by  our  new  method  described  in  reference  (a) . 

Acquisition  probabilities  are  computed  and  compared  with  those 
obtained  in  reference  (a)  with  dog-leg  unmasking. 


This  research  was  supported  by  the  Office  of  Naval  Research 
under  Contract  No.  N00014-70-C-0322. 
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Summary  and  Conclusions 

The  main  conclusion  reached  from  the  study  is  that  our  new  TMA 
and  the  associated  method  of  computing  lead  angles  will  counteract 
the  adverse  effects  of  torpedo  radiated  noise,  provided  that  the 
sonar  tracks  the  centroid  of  the  received  signals,  and  the  signals 
are  merged  sufficiently  (both  signals  in  the  main  tracking  lobe) 
when  post- launch  tracking  begins.  By  simply  tracking  the  centroid 
bearing,  accurate  torpedo  tracks  for  interception  are  generated 
in  the  corrected-intercept  control  mode,  even  when  the  target 
maneuvers  radically  and  the  torpedo  noise  intensity  at  the  source 
exceeds  the  target  sound  intensity  at  the  source.  Acquisition 
probabilities  are  comparable  with  those  obtained  in  reference  (a) 
with  dog- leg  unmasking. 

The  good  results  that  are  obtained  with  the  corrected-intercept 
control  mode  from  bearing  data  obtained  in  centroid  tracking  will 
not  be  obtained  with  other  known  methods  of  target  motion  analysis, 
such  as  CHURN,  and  the  associated  methods  of  computing  fire 
control  orders.  The  favorable  outcome  stems  from  the  motion  of 
the  centroid  bearing  line,  the  effect  of  this  motion  on  the 
apparent  bearing  rate,  and  the  strong  use  of  current  bearing 
rate  to  compute  the  lead  angle  in  our  solution. 

Centroid  Tracking 

We  assume,  as  before,  that  the  target  is  completely  masked  ar.u 
tracking  is  impossible,  when  the  torpedo  is  very  close  to  the 
tracking  submarine.  If  desired,  we  can  use  the  square  spreading 
law  to  express  the  range  at  which  complete  masking  terminates  in 
terms  of  the  range  of  the  target,  the  difference  in  signal 
strengths  at  the  source,  the  difference  in  sonar  sensitivities 
in  the  two  directions,  and  the  difference  in  signal  strengths  at 
the  receiver  that  corresponds  to  the  end  of  complete  masking. 

Thus,  let 
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AS  =  difference  in  received  signal  strengths  (target  minus  torpedo) 
AH  =  corresponding  difference  in  signal  strengths  at  the  source 
AG  =  corresponding  difference  in  directional  sensitivities 
R  =  range  of  the  target  from  the  tracking  submarine 
R^,  =  range  of  the  torpedo  from  the  tracking  submarine 


Then 


AS  =  AG  +  AF  -  20  log1Q  (R/R^,) , 


(1) 


from  which  we  obtain 

_  _  0.115  (AS  -  AG  -  AH) 

RT  =  R  e 


(2) 


Here  AS,  AG,  and  AH  are  in  decibels,  and  AG  is  positive  when  the 
target  is  in  the  major  lobe  and  the  torpedo  is  in  a  side  lobe. 

In  equation  (1)  we  have  omitted  the  attenuation  term  because  it 
is  small  relative  to  other  terms  at  the  ranges  of  interest. 

To  use  equation  (2)  to  find  the  range  at  which  complete  masking 
ends  we  need  to  specify  a  value  for  AS  that  corresponds  to  this 
condition.  The  value,  AS  =  0,  usually  is  considered  to  be  the 
value  at  which  the  target  signal  ceases  to  be  merged  with  the 
torpedo  signal,  so  that  the  two  signals  might  be  tracked 
separately.  However,  it  is  difficult  to  estimate  the  value  of 
AS  that  represents  the  passage  from  no  tracking  to  "some" 
tracking  of  the  target. 

We  can  avoid  this  difficulty  by  assuming  that  the  sonar  operator 
always  tracks  a  composite  signal  at  a  bearing  that  is  a  weighted 
average  of  the  two  separate  bearings.  Let 


B  =  bearing  of  the  target 
Bt  =  bearing  of  the  torpedo 
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The  apparent  bearing  of  the  composite  signal  is  the  weighted 
average, 


B*  -  (B  +  oB^)  /  (1  +  c), 
where 

p  =  10-iS/1° 


(3) 


£ 


Substituting  equation  (1)  into  equation  (-5)  ve  obtain 


P  -  (R/R /  e-°-23  (-G  +  iH> 


(5) 


If  the  torpedo  and  target  submarine  are  in  the  major  lobe;,  ve 
usually  assume  that  AG  =  0.  This  assumption  is  acceptable  when 
the  sensitivity  is  nearly  constant  over  a  large  fraction  of  the 
lobe.  If  the  major  lobe  is  sharply  peaked,  better  estimates  of 
AG  and  p  can  be  obtained  from  the  sensitivity  function.  Thus, 
let  G(b)  be  the  sensitivity  at  angle  b  from  the  axis.  Then 


AG  =  G(B  -  B*)  -  G(Bt  -  B*) 

Using  equation  (3),  AG  becomes 

AG  =  GfoAB/U  +  p))  -  G(-AB/(1  +  o» ,  (6) 

where 

AB  =  B  -  Bt  (?) 

We  then  find  the  values  of  o  and  AG  from  equations  (5)  and  (6)  by 
iteration.  Start  with  AG  =  0,  find  o  from  (5)  and  substitute  in 
(6)  to  obtain  a  new  value  of  AG^  repeat  until  no  change  Is  obtained. 

From  equations  (3)  and  (A)  it  is  seen  that  o  =  1  and  the  bearing  angle 
B*  is  the  mean  of  the  two  bearing  angles  when  AS  =  0.  Also,  it 
is  seen  from  (6)  that  AG  =  0  when  p  =  1,  If  the  function  G(b)  is 
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symmetrical  about  the  axis,  hence,  if  we  put  .IS  —  0  when  htch 
target  and  torpedo  are  in  the  major  lobe,  the  estimate  of  e  will 
be  accurate  in  the  vicinity  of  the  cross-over?  it  -will  be  too 
large  when  c  <  1  and  too  small  when  c  >  1,  That  is,  the  apparent 
bearing  line  will  be  closer  to  the  mean  bearing  line  when  !u5  is 
set  equal  to  zero  than  would  be  the  case  when  13  is  estimated  from 
the  sensitivity  function. 

For  example,  consider  the  sensitivity  function 

G(b)  =  G  -F  G  cos  (~b/2B’ ) ,  b  <B' 

1  o 

=  Gx  ,  b  ^ 

where  £'  is  the  half-angle  of  the  major  lobe,  S.  is  the 

sensitivitv  in  the  minor  lobes,  and  G  is  the  tain  in  ssnsitivitv 
"  o 

on  the  axis  of  the  main  lobe  over  that  in  the  minor  lobes,  it 
both  the  target  and  the  torpedo  are  in  the  major  lobe, 

dG  =  G  f cos  (r£/2(l  +  p-1))-  cos  (r£/2&  05 

o 

where 

£  =  >b:  /  S1  0.2)5 

The  expression  (9)  for  1G  is  obtained  directly  from  equations  05 
and  (S).  It  applies  when  both  arguments  are  less  than  or  equal 
to  «/2,  which  is  equivalent  to  the  condition 

S  <_  1  +  min  (c,c  ill) 

It  is  evident  that  1G>0  when  oO,  and  1G<0  when  eO. 

We  estimate  p  by  iteration  from  equations  (55  and  (-95.  The  first 
approximation  is 
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(12) 


/D/D  x-  -0.23AH 
^  =  (R/Rj.)  e 

with  AG  =  0.  The  next  estimate  of  AG  is 

A^G  =  AG  in  (9)  with  p=o^  (13) 

and  the  corresponding  approximation  for  p  is 

,  .  ,x  e-°-23i1G  (U) 

If  condition  (11)  'is  not  satisfied,  the  sonar  either  is  tracking 

the  target  or  the  torpedo,  depending  on  the  signal 

strengths,  the  skill  of  the  operator,  and  other  circumstances. 

In  general,  if  p^«l,  it  is  the  target;  and  if  P^>>1,  it  is  the 

torpedo.  For  close  to  1,  either  one  could  be  tracked,  depending 

on  which  one  happened  to  be  in  the  major  lobe  when  the  separation 

between  the  bearing  lines  became  large  enough  that  condition  (11) 

is  not  satisfied. 

\ 

Condition  (11)  certainly  is  satisfied,  if  8<1;  and  is  not  satisfied, 
if  iv2.  For  1<3<2,  the  condition  is  satisfied  for  some  values  of 
c  and  is  not  satisfied  for  others.  We  must  replace  the  cosine 
term  for  which  the  inequality  is  not  satisfied  by  0.  Thus,  if  we 


Y  -  min  [1,  i$/(l  +  o)],  y  =  min  [1,  pg/(l  +  p)]. 


the  equation 

AG  =  Gq  [cos  (~y^/2)  -cos  (ny/2) J 


applies  for  S^2.  In  making  the  approximation,  we  put  p  =  p^  in 

1  ] 

(15),  compute  A^G  by  putting  the  corresponding  values  of  y  and  y 
in  (16),  and  then  substitute  in  (14)  to  estimate  p. 
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If  S>2»  it  should  be  possible  to  detect  the  separate  signals.  If 
there  is  any  doubt  as  to  which  one  is  the  target,  the  doubt  can 
be  resolved  by  dead-reckoning  the  torpedo.  Hence,  we  assume  that 
the  major  lobe  will  be  put  on  the  target,  and  that  consequently 
AG  =  Gq  for  all  values  of  o  •  hen  £>2.  This  assumption  might 
be  questionable  when  p^>>1,  since  then  the  received  target  signal 
with  no  gain  is  much  weaker  than  the  corresponding  received 
torpedo  signal.  A  possibility  exists  that  under  these  conditions  - 
and  under  some  other  conditions,  as  well  -  the  operator  may  track 
on  the  wrong  null  in  the  returned  target  signal,  which  would 
introduce  a  large  bias.  Our  acquisition  models  do  not  include 
this  possibility. 

Bearing  Rates 

The  motion  of  the  centroid  bearing  line  depends  on  the  motions 
of  the  target  bearing  line  and  the  torpedo  bearing  line,  and  on 
the  rate  of  change  of  p.  From  equation  (3)  we  obtain 

B*  =  (B  +  p  Bx)  /  (1  +  p)  -  \B  p/(l  +  p)2  (17) 

The  first  term  on  the  right-hand  side  of  equation  (17)  is  the 
weighted  mean  bearing  rate,  while  the  second  term  Is  produced 
by  the  changing  weight.  It  was  expected  that  the  first  term 
would  be  the  dominant  term,  in  general.  However,  In  a  sample 
computation  the  magnitude  of  the  second  term  was  larger  than 
that  of  the  first  term  at  the  start  of  the  post-launch  tracking, 
when  |AB[  was  large.  This  term  has  the  same  sign  as  AB,  since 
p  is  negative. 

We  can  write  p  in  the  form 

P  =  gh  (R/Rj)2  (IS) 
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where 


-.23AG  .  --232H 

g  =  e  ,  h  -  e 


(19) 


and  AG  is  the  difference  in  sensitivity  gain  and  13  is  the 
difference  in  signal  strengths  at  the  source.  The  factor  h  is 
constant  and  the  factor  g  decreases  as  t  increases,  since  IG 
increases.  Also,  R/R decreases  as  t  increases.  From  (18)  ve 
have 


c  -  h  g  (R/Rt)2  +  2  gh  (R/rJ)  (R_R  -  ai.)  (20) 

Both  terms  on  the  right-hand  side  are  negative,  since  g  <  0, 

•  0 

R  <  R^.,  and  R^  <  R  for  the  control  interval.  The  g  tent  may  fee 
significant  at  the  start  of  tracking,  as  in  our  example;  it  is 
difficult  to  express  in  term  of  the  rates  of  change  of  the 
bearings  and  ranges,  since. g  depends  on  c  in  a  transcendental 
form.  At  any  rate,  o  usually  decreases  fast  enough  to  give  3* 

.in  equation  (17)  a  large  pseudo  component  toward  the  bearing  3 
of  the  target. 

The  Role  of  the  TMA  and  Correction  Computations 

The  effect  of  the  large  component  of  pseudo  bearing  rate,  in  the 
direction  from  the  torpedo  bearing  toward  the  target  bearing,  will 
depend  very  strongly  on  the  target  notion  analysis.  For  a  method, 
such  as  CHURN,  in  which  the  target  parameters  are  obtained  from 
a  simultaneous  solution  based  on  all  the  observed  bearings,  the 
effect  of  a  small  number  of  additional  bearing  observations  after 
launch  will  have  a  small  effect  on  the  solution.  The  effect 
will  be  large,  and  difficult  to  predict,  when  the  number  of  post¬ 
launch  observations  is  comparable  to  the  nucber  made  before 
launch;  it  will  depend  on  the  target  motion. 
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For  example,  consider  the  target  motion  in  the  check  calculation, 
run  3/case  3  of  reference  (a) :  the  target  runs  on  course 
90  degrees  at  5-6  yds/sec.  for  400  seconds  and  then  turns 
120  deg;rees  to  course" -30  degrees.  The  bearings  of  the  target 
and  torpedo,  relative  to  the  position  of  own  submarine  at -the 
start  of  tracking  (that  is,  subtracting  off_the  effect  of  own 
submarine  motion  on  the  bearings) ,  are  shovm  in  Figure  1.  Also 

snown  are  the  centroid  bearing  'and  the  torpedo  bearing  after 

2  ■ 

post-launch  tracking  begins,  assuming  ho  corrections  are  made. 

•k  , 

The  tracks  of  the  torpedo  bearing  and  the  centroid  bearing  B 
will  depend  on  the  corrections  that  are  made  to  the  torpedo  gyro 
course,  and  the  corrections  will  depend  on  the,  TMA.  If  the 

*  <7 

CHURN  IMA  is  used,  there- will  be  little  change  in  B^.  and  B  for 
the  first  few  Opportunities,  since  wc  have  assumed  that  corrections 
will  be  made  at  intervals  of  20  seconds,1  which  /s  only  5  percent- 
of  the  pre-launch  tracking  interval.  When  the  post-launch  tracking 
interval  becomes  appreciable  -  say,  25  percent  or  more  -  it  is 
difficult  to  anticipate  what  the  CHURN  solution  will  do  to  the 
torpedo  course,  without  making  a  detailed  computation.  Certainly, 
when  the  post- launch  tracking  interval  is  at  least  50  percent 
that  of  the  pre-launch  interval,  CHURN  will  have  a  difficult 
problem  in  trying  to  fit  a  single  tangent  curve  to  these  bedring 
observations .  In  fact,  the  problem  would  be  difficult  even  if  the 

*  V 

target  did  not  maneuver. 

A  similar  result  is  obtained  for  most  methods  of  target  motion 
analysis,  including  the  equal-segment  method,  and  the  partially- 
mechanized  version  of  it  known  as  MATE.  A  few  observations' in 
post— launch  tracking  will  have  a  small  effect.  The  effect  of  a 
large  number  of  observations  in  post-launch  tracking  will  depend 
on  the  relative  weights  that  are  given  to  the  post-launch  data 
and  the  pre- launch  data. 

The  motion  of  the  torpedo  when  corrections  are  computed  from  our 
IMA  is  quite  different.  At  the  first  opportunity  a  large 
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correction  is  applied  in  the  direction  toward  the  target. 

Additional  corrections  in  this  direction  are  applied  at  later 
opportunities,  until  the  torpedo  bearing  is  close  to  the  target 
bearing  and  |£b|  is  snail.  In  the  absence  of  danping,  the 
torpedo  usually  swings  to  the  opposite  side  of  the  target  bearing 
line  before  corrections  are  applied  to  drive  it  back.  Thus, 
oscillations  occur,  unless  sone  fom  of  danping  is  applied. 

Tteo  sinple  foms  of  danping  were  found  to  be  effective.  The 
first  one  often  is  called  proportional  navigation  in  guided- 
.  nissile  applications .  The  second  one  is  a  sinple  linitation  on 
the  nagnitude  of  the  change  in  course  angle  that  can  be  made  at 
any  correction.  Kith  either  forn  of  danping  the  oscillations 
can  be  controlled  -  even  eliminated  entirely,  if  desired  -  to 
obtain  accurate  interception  courses. 

The  reason  for  the  fast  reaction  to  the  apparent  target  notion 

in  the  post-launch  phase  is  the  emphasis  put  on  the  current 

bearing  rate  in  our  TMA  and  corresponding  solution. for  the  lead 

angle.  We  use  only  the  post-launch  data  in  computing  the  bearing 

rate,  and  we  use  exponential  weighting  to  keep  it  current.  Ke 

* 

then  compute  the  component  u^  of  target  velocity  normal  to  the 
sight  line,  and  the  corresponding  lead  angle,  using  the 
estimate  of  range  made  at  launch  and  assuming  that  the  range 
component  of  velocity  is  zero. 

•  #  / 

The  large  component  of  B*  obtained  from  the  p  term  in  equation 

(17)  may  yield  unrealistic  values  for  for  a  few  corrections. 

(In  our  example  we  obtained  a  value  of  approximately  -  30  yds./ 

sec.  for  two  corrections.)-  To  avoid  over-correcting  and  to  keep 

the  computer  happy,*  we  restrict  |u^|  to  be  no  more  than  some 

Unless  ju^|  is  less  than  the  torpedo  speed,  the  formula  used  in 
computing  the  lead  angle  0  may  produce  a  value  of  sin  6  that  exceeds 

( •  L  '  L 

1.0  in  absolute  value,  which  causes  the  computer  to  complain  and 
sometimes  to  react  in  a  spiteful  way. 
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reasonable  limit,  such  as  15  yards/sec.,  that  is  less  than  the 
torpedo  speed. 

An  attempt  was  made  to  study  the  motion  of  the  torpedo  by  writing 
the  differential  equations  of  motion,  under  the  assumption  that 
corrections  are  applied  continuously.  The  equations  are 
difficult  to  solve  analytically.  A  numerical  solution  could  have 
been  obtained  by  a  number  of  methods,  but  it  would  have  been  a 
less  accurate  description  of  the  true  motion  than  that  which  we 
obtain  from  the  deterministic  computer  simulation  that  is 
embedded  in  our  analytic  simulation  model. 

Centroid  Capture  by  the  Torpedo 

An  important  effect  of  the  p  term  in  B*,  when  B*  is  used  in  our 
TMA,  is  to  avoid  having  the  centroid  bearing  lock  onto  the 
torpedo  bearing.  The  biggest  danger  occurs  at  the  start  of 
post-launch  tracking. 

If  |B*  -  B|  decreases  at  the  star.t  of  post-launch  tracking,  centroid 
capture  by  the  torpedo  wilT  hot  occur  immediately  and  is  not 
likely  to  occur  later.  Now,  | B*  —  B |  will  decrease,  if  B*  -  B 
has  the  same  sign  as  B  -  B*,  which  is  the  same  as  the  sign  of 
B  -  B^.  Hence,  to  avoid  capture  we  want  (§*  -  6)  /  (B  -  B^)  >  0. 
From  equation  (17)  we  obtain 

B*  -  B  _  -  P  (  _ I 

b-bt  1  +  P  \B_Bt 

Hence,  we  have  the  condition 

'  {nh 

or 

~  log  ( (1  +  1/  p)  /  |  AB  | ]  >  0, 
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Arthur  I)  little  Inc 


which  is  satisfied  if  and  only  if 
(1  +  1/p)  /  | ABj 

is  an  increasing  function  of  t.  Since  the  numerator  is  an 
increasing  function  of  t,  the  only  danger  of  capture  is  when 
|ab)  is  increasing,  and  at  a  fractional  rate  that  is  greater 
than  that  of  (1  +  1/p). 


The  Simulation  Model 

The  simulation  model  for  the  corrected-intercept  mode  that  is 
described  in  reference  (a)  was  revised  to  remove  the  dog-leg 
unmasking  and  to  insert  the  simulation  of  centroid  tracking.  The 
required  changes  in  the  computation  of  the  acquisition  probability 
also  were  made.  The  revised  computational  procedure  is  outlined 
in  the  Appendix,  with  only  the  revised  sections  written  in  detail. 


The  rev'sed  procedure  for  centroid  tracking  was  programmed. 
Acquisition  probabilities  were  computed  for  the  same  runs  used  in 
reference  (a).  Some  of  the  new  parameters  were  varied  to 
determine  their  effects,  but  most  of  them  were  assigned  one 
arbitrary  value.  The  results  are  given  below. 


Comparison  Conditions:  Runs  and  Parameter  Values 

The  five  run  types  are  maneuvers  from  an  initial  course  of 
90  degrees  and  target  speed  of  10  knots,  starting  at  a  range  of 
10,000  yards.  The  maneuvers  are  as  follows: 


Run  Maneuver _ 

1  60  degrees  turn  away 

2  60  degrees  turn  towards 

3  120  degrees  turn  away 

4  120  degrees  turn  towards 

5  Decelerate  from  10  knots  to  4  knots  and 

accelerate  back  to  10  knots 
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Own  submarine  runs  at  5  knots  on  alternate  courses  of  90  degrees 
and  -90  degrees  for  tracking  legs  of  100  seconds,  with  50  seconds 


allowed  for 

times: 

Case 

the  turns. 

Time 

The  maneuvers  are  started  at  the  following 

Position  of  Own  Submarine 

1 

200  secs. 

Middle  of  second  tracking  leg 

2 

300  secs. 

Start  of  third  tracking  leg 

3 

400  secs. 

End  of  third  tracking  leg 

4 

500  secs. 

Middle  of  fourth  tracking  leg,  if  made 

5 

600  secs. 

Start  of  fifth  tracking  leg,  if  made 

The  tracking 

and  fire- 

control  procedure  described  in  Reference  (a) 

uses  tests  at  the  end  of  the  first  three  tracking  legs  to 
determine  whether  or  not  to  accept  the  range  solution  that  has 
been  generated  and  launch  the  torpedo,  or  to  track  for  an 
additional  leg.  The  test  consists  essentially  of  a  comparison 
of  our  new  estimate  r  based  on  three  legs  with  the  Ekelund 
estimate  r  based  on  the  last  two  legs.  It  is 

r.  <r<r  ,  r  .  <r<r 
mxn  max  mm  max 

a  <  r/r  <  1/a  ,  0  <  a  <1 
r  r  r 

We  test  to  see  if  both  estimates  are  in  a  reasonable  range 

interval  and  are  close  together.  If  all  inequalities  are  satisfied, 

we  accept  the  estimate  r.  If  any  inequality  is  not  satisfied, 

we  track  for  an  additional  leg  and  try  again.  Simulation  runs, 

reported  in  Reference  (a),  show  that  the  test  usually  rejects 

a  poor  solution  and  usually  accepts  a  good  solution  when  a^  =  0.7, 

r  =  1000  yards,  r  =  25,000  yards.  The  values  of  r  .  and 
min  max  min 

r  are  not  critical, 
max 

If  the  solution  at  the  end  of  three  legs  is  accepted,  as  usually 
occurs,  the  maneuver  is  made  during  tracking  in  cases  1  and  2,  at 
the  end  of  tracking  in  case  3,  and  after  tracking  has  terminated 
in  cases  4  and  5. 
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The  parameter  values  listed  in  Table  3.1  of  reference  (a)  were 

retained,  except  that  r  .  =  2  d  was  used  to  avoid  the  slight 

r  min  g  & 

possibility  of  accepting  a  range  estimate  less  than  d  ;  this 

6 

change  had  no  effect  on  the  comparisons. 


Some  of  the  parameters  listed  in  Table  3.1  of  reference  (a)  were 
eliminated  when  the  dog  leg  was  eliminated,  and  other  parameters 
were  introduced  in  the  centroid  tracking.  Additional  parameters, 
mostly  new,  and  the  values  assumed  for  them  in  the  initial 
comparison,  are  listed  in  Table  1.  After  making  computations 
for  all  runs  and  cases  with  the  parameter  values  listed  in 
Table  1,  additional  computations  were  made  to  explore  the  effects 
of  some  of  the  parameters,  such  as  AH. 


TABLE  1.  ADDITIONAL  PARAMETER  VALUES 


Parameter 

Units 

Value 

Definition 

fl 

fraction 

0.4 

Start  post-launch  tracking: 
fraction  of  estimated  range 

f2 

fraction 

0.5 

Enable:  fraction  of  estimated 
range  (new  definition) 

AH 

decibels 

5.0 

Target  sound  above  torpedo  noise 

G 

0 

decibels 

14.0 

Gain  at  center  of  major  tracking 
lobe  over  side  lobes 

B1 

radians 

0.3 

Half  angle  of  major  tracking  lobe 

B1 

ratio 

2.0 

Discrimination  ratio  of  separation 
angle  to  the  angle  B^ 

V 

fraction 

0.4* 

Proportion  of  correction  in 
Proportional  Navigation 

A0 

max 

radians 

0.2* 

Limit  on  corrections 

max  |uj 

yds/sec. 

15.0 

Limit  on  normal  velocity  component 

not  used  simultaneously;  the  set  (f^  =  1.0,  &®maxi=  0*2)  and 

the  set  (f.  =  0.4,  A6  =  1.0)  were  used. 

4  max 
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Another  fraction  f^,  with  a  value  of  0.9,  had  been  inserted  to  cut 
off  the  control  corrections,  to  avoid  wild  oscillations  near 
expected  passage.  However,  after  damping  (proportional  navigation 
or  limited  corrections)  had  been  applied,  it  was  found  that 
control  cutoff  was  not  needed  and  sometimes  stopped  corrections 
too  soon.  Hence,  it  was  removed,  in  effect,  by  putting  f^  =  20. 

Before  choosing  values  for  the  damping  parameters,  f^  and  , 

a  partial  sensitivity  analysis  was  made  for  run  3/case  3.  With 

A0  =  1.0,  f,  was  varied  from  0.2  to  0.6  inclusive;  with  f.  =  1.0, 
max  4  *  4 

A0  was  varied  from  0.1  to  0.3  inclusive.  The  large  value  of 
max 

the  parameter  that  was  not  varied  had  the  effect  of  removing  the 
corresponding  form  of  damping.  For  both  forms  of  damping  the 
acquisition  probability  remained  almost  constant  when  the  damping 
parameter  was  varied,  indicating  that  the  chosen  interval  was 
on  the  plateau  of  the  curve,  since  very  low  acquisition 
probabilities  were  obtained  with  no  damping  and  with  excessive 
damping.  From  these  results  the  following  two  sets  of  damping 
parameters  were  chosen: 

Proportional  Navigation:  f^  =  0.4,  A3^ay  =  1.0 

Limited  Correction  :  f .  =  1.0,  A3  -  0.2 

4  max 

In  each  case  the  correction  applied  to  the  gyro  course  is 

+  f,  min  (|a9|,  AS  ),  if  A3>0 
4  1  1  max  — 

-  f.  min  (I  AS  I,  AS  ),  If  A3<0 
4  max 

where  A0  is  the  correction  that  is  computed  from  the  IMA.  If 

proportional  navigation  alone  is  used,  the  correction  is  simply 

f ,  AO . 

4 
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Results 


Some  results  are  shown  in  Tables  2  and  3.  Acquisition  probabilities 
for  all  runs  and  cases  are  shown  in  Table  2,  for  rbe  parameter 
values  listed  in  Table  1,  including  the  two  forms  of  damning-. 

Also  shown  for  comparison  purposes  are  the  acquisition 
probabilities  from  Table  3. A  of  reference  (a)  when  no  -control  is 
exercised  and  when  a  dog  leg  is  used  to  unmask. 

Acquisition  probabilities  are  shown  in  Table  3-,  for  runs  3  and  4 
onlv ,  with  AH  equal  to  5,  0,  and  -5  decibels,  to  show  the  effect 
of  changes  in  the  target  sound  intensity  and  the  torpedo  noise. 

The  value  of  AH  is  the  sound  intensity  of  the  target  above  the 
torpedo  noise. 
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TABLE  2.  ACQUISITION  PROBABILITIES-CORRECTED  INTERCEPT  MODE 


Run 

Case 

ho  Control 

Centroid  Tracking: 
Dog  Leg  40%  Correction  12°  max 

AH  =  5 

;  Correction 

1 

1 

.34 

.46 

.45 

.43 

2 

.49 

.56 

.52 

.50 

3 

.47 

.54 

.57 

.55 

4 

.50 

.57 

.61 

.60 

5 

.54 

.59 

.64 

.64 

r 

1 

.77 

.79 

.82(a) 

.81(a) 

2 

.75 

.79 

.81 

.80 

3 

.69 

.75 

.79 

.79 

4 

.67 

.74 

.77 

.77 

5 

.65 

.72 

.75 

.75 

5 

1 

.32 

.42 

.41 

.41 

? 

.34 

.46 

.45 

.45 

3 

.33 

.51 

.56 

.55 

4 

.42 

.50 

.62 

.60 

5 

.51 

.53 

.66 

.66 

4 

1 

.80 

.84 

.86 

.86 

2 

.78 

.84 

.86 

.86 

3 

.58 

.74 

.73 

.73 

4 

.63 

.71 

.72 

.75 

5 

.65 

.66 

.73 

.74 

5 

1 

.59 

.67 

.67 

.68 

2 

.58 

.61 

.68 

.67 

3 

.59 

.68 

.68 

.69 

4- 

.59 

.68 

.69 

.69 

5 

.59 

.68 

.69 

.69 

Soter 

(a)  These  acquisition  probabilities  are  obtained  when 

the 

torpedo  is  enabled  at  the'  time  ,at  which  the  laminar  point 
has  reached  the  50  percent  point  (f ^  =0.5)  of  the  estimated 
range.  If  enable  is  delayed  until  f2  =  0.7,  the  laminar 
point  has  almost  certainly  passed  the  target  and  the 
acquisition  probabilities  are  very  small.  See  the 
section.  Discussion  of  Results,  for  an  explanation. 
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TABLE  3.  ACQUISITION  PROBABILITIES:  EFFECTS  OF  TORPEDO 


AND  TARGET  NOISE 


40%  Correction  12  max.  Correction 


Run 

Case 

AH(a)-5 

0 

-5 

5 

0 

3 

1 

.41 

.50 

.50 

.41 

.50 

.49 

2 

.45 

.49 

.48 

.45 

.48 

.37 

3 

.56 

.55 

.49 

.55 

.55 

.55 

4 

.62 

.61 

.60 

.60 

.60 

.60 

5 

.66 

.64 

.65 

.66 

.62 

.63 

4 

1 

.86 

.86 

.84 

.86 

.80 

.64 

2 

.86 

.85 

.84 

.86 

.83; 

.67 

3 

.73 

.80 

.73 

.73 

.78' 

.79 

4 

.72 

.76 

.79 

.75 

.70 

.75 

5 

.73 

.72 

.76 

.74 

.70 

.64 

Note: 

(a) 

AH  =  Target 

sound 

intensity 

at  source 

minus  Torpedo 

noise  at  source 
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Discussion  of  Results 


From  Table  2  it  is  seen  that  when  the  additional  parameters  have 
the  values  listed  in  Table  1  the  acquisition  probabilities  with 
centroid  tracking  of  the  degraded  signal  are  comparable  with 
those  obtained  by  unmasking.  In  some  cases,  such  as  cases  4  and 
5  of  run  3,  they  are  significantly  higher.  These  results  suggest 
that  centroid  tracking  probably  is  significantly  better  than 
dog-leg  unmasking  when  the  target  turns  away  after  being  alerted 
by  the  launching  noises,  and  nay  be  socewhat  better  against  most 
maneuvers  at  longer  ranges.  These  questions  can  be  explored  by 
more  extensive  computations. 


It  also  appears  from  Table  2  that  there  is  little  to  choose 
between  the  two  forms  of  damping  with  the  values  chosen  for  the 
damping  parameters.  Results  in  Table  3  suggest  that  proportional 
navigation  is  slightly  better  than  limited  correction:  However, 
the  values  chosen  for  the  damping  parameters  are  not  necessarily 
optimal  for  a  particular  run  and  case,  or  for  a  mixture  of  runs 
and  cases.  Again,  more  extensive  computations  are  needed  to 
answer  these  questions. 


The  enable  fraction  f ^  was  first  chosen  to  be  0.7  to  delay  enable 
until  the  laminar  point  is  70  percent  of  the  distance  to  the 
expected  interception  point.  For  case  1  of  run  2  the  acquisition 
probability  was  0.01  for  both  forms  of  damping.  Examination  of 
the  details  revealed  that  enable  had  been  delayed  too  long  and 
the  laminar  point  had  passed  the  target  (with  high  probability) 
when  enable  occurred.  For  this  type  of  maneuver  and  the  assumed 
path  of  own  submarine  during  tracking,  both  range  estimates) 
r  and  r  are  too  large  when  three  tracking  legs  are  used,  but  are 


close  enough  together  to  be  accepted  by  our  test-  The  70  percent 
value  of  the  accepted  range  estimate  exceeds  the  true  range-  to 
passage  by  the  laminar  point,  particularly  since  the  target  path 
is  such  as'  to  decrease  the  range.  In  case  2  of  run  2  enable  occurs 
barely  in  tine. 
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r  i 


r  » 
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The  late  enable  time  that  occurred  in  case  1  of  run  *2  does  not 
occur  in  case  1  of  run  4,  even  though  the  target  again  turns 
toward  own  submarine  at  the  execution  time  of  200  seconds.  The 
120  degrees  change  in  course,  rather  than  the  60  degrees  in  run  2, 
produces  underestimates  of  range  by  both  methods.  The  errors  are 
large,  but  the  tvo  estimates  again  are  close  enough  to  be  accepted. 
An  understanding-  of  this  paradox  can  be  obtained  by  a  detailed 

i 

examination  of  bearing  rates  on  the  tracking  legs  and  their 
effects  on  the  range  estimates. 

v  ,  <  l 

To  avoid  late  enable,  such  as  that  which  occurred/ 'in  case  1  of 

// 

•  run  2,  we  reduced  f ^  to  0.5  fqr  these  computations.  It  is  .  ' 

questionable  whether  the  avoidance  of  an  occasional  late  enable 
more  than  offsets  the  delay  in  passage. that  would  occur  when  the 
torpedo  runout  speed  exceeds  the  search  speed.  The  question  is 
/^di-fficult  to  answer,  even  with  extensive  computations,  since 
it  involves  distributions  of  the  engagement  range  and  the  target 
course  and  speed.  . 

A  natural  question-  is  the  extent  to  which  the  results  depend  on 

the  assumed  values  for  the  parameters  that  control  the  centroid 

tracking;  particularly  the  difference  AH  in  sound  output  from 

the  target  and  torpedo.  The  comparison  in  Table  3  shows  that 

»  * 
there  is  little  change,  and  not  always  downward,  in  the 

acquisition  probability  on  runs  3  and  4  for  either  form  of 

damping  when  AH  changes  from  5  db  to  0  db.  With  proportional 

navigation/there  are  no  significant  changes  when  AH  is  dropped 

to  -5  db ^  With  the  12  degrees  limitation  on  the  correction, 

decreases  occur  in  a  few  cases. 
l 

i  > 

In  some  cases,  such  as  case  1  of  run  3,  the  torpedo  captures 
the  tracking  beam.  For  this  case  the  range  estimates  r  and  r 
at  the  end  of  three  legs  are  only  about  40  percent  of  the  true 
range,  and  are  close  enough  together  to  be  accepted.  With  the 
resulting  small  range  estimate  post-launch  tracking  starts  soon 
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after  launch,  and  the  value  of  the  weighting  factor  p  is  very 
large,  even  for  AH  =  5  db.  Since  |AB]  is  increasing  and  at  a 
fractional  rate  that  exceeds  that  of  (1  -r  1/p),  the  torpedo 
captures  the  tracking  bean  at  the  start  of  post-launch  tracking. 
However,  the  large  value  of  p  eventually  contributes  a  sufficient 
component  to  the  bearing  rate  in  the  direction  of  the  target  to 
generate  corrections  to  the  gyro  course  that  break  the  capture 
condition. 

The  conditions,  if  any,  under  which  capture  persists,  after  it 
occurs,  are  difficult  to  determine,  since  the  condition  depends 
in  a  complex  way  on  the  motions  of  the  target  and  torpedo  relative 
to  own  submarine,  the  range  estimate,  the  ranges  to  the  torpedo 
and  target,  and  many  control  parameters.  The  critical  question 
is  whether  the  early  course  "corrections"  that  are  generated 
move  the  torpedo  bearing  closer  to  the  target  bearing  or  farther 
away.  The  question  of  the  conditions  under  which  torpedo  capture 
of  the  tracking  bean  persists  is  an  important  question  in  the  use 
of  corrected-intercept  control  with  our  TMA  against  targets  that 
are  quieter  than  the  torpedo.  If  the  target  is  quieter  than  the 
torpedo,  the  torpedo  usually  captures  the  tracking  bean  at  the 

4  . 

start  of  post-launch  tracking.  The  use  of  our  TMA  apparently 
permits  capture  to  be  broken  before  the  laminar  point  passes  the 
target. 

When  the  separation  angle  |ab{  becomes  large  enough  the  sonar 

operator  can  detect  the  separate  signals  and  presumably  shift 

the  tracking  beam  to  the  target.  An  allowance  is  made  for  this 

shift  in  our  model:  if  6  >  b\  the  gain  G  in  sensitivity  shifts 

to  +  G  from  the  value  near  -  G  it  had  when  the  torpedo  was 
o  1  o 

tracked.  Here,  6  =  |CB|  /B  is  the  half  angle  of  the  tracking 
lobe.  We  assume  that  the  operator  can  discriminate  between  the 
two  signals  when  the  separation  angle  is  6^  times  the  half  angle. 
In  our  computations  we  used  B^  =  0.3  radians  and  6^  =  2.0,  which 
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is  equivalent  to  a  separation  requirement  of  34  degrees.  A 
smaller  separation  requirement  could  permit  a  shift  to  be  made 
sooner.  The  effect  on  the  acquisition  probability  depends  on 
many  factors.  A  shift  in  the  value  of  to  0.15  radians  yielded 
lower  acquisition  probabilities  for  case  1  of  run  3. 


Ref erer 

(a)  ADL  Report  XUKRES  ••'12 ,  "Control  Modes  and  Acquisition 
Probabilities  for  Torpedo  MX  48  (U) ,  "Contract  No. 
N00140-68-C-0278,  January  1970,  Final  Report  Confidential, 
Technical  Appendices  (bound  separately)  Unclassified. 
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APPENDIX  A 


COMPUTATIONAL  PROCEDURE  FOR  THE  C-I  MODE  WITH  CENTROID  TRACKING 


The  computational  procedure  for  the  corrected- intercept  mode  when 
a  dog  leg  is  used  to  unmask  is  given  in  Appendix  E  of  Reference  (a) . 
We  list  below  only  the  changes  that  are  needed  to  replace  the 
unmasking  maneuver  by  centroid  tracking. 


1.  TMA  for  Torpedo  Launch 


2.  Initial  Course  and  Launch 

These  two  sections  are  not  affected  by  the  masking  assumptions, 
and  are  equivalent  to  section  E.2  of  Appendix  E  of  Reference  (a). 


3.  Post-Launch  Tracking 

We  start  with  the  coordinates  (x  y  )  of  the  pseudo  origin, 

io,  lo  ^ /  \ 

the  initial  gyro  angle  0  ,  the  estimate  rv 1  of  range,  and  the 
/  \  o 

number  Nv  of  time  steps  in  the  pre-launch  tracking  interval. 


f^  =  factor  for  initial  leg  (input) 

f^  =  factor  for  the  enabling  step  (input) 

=  torpedo  speed  during  runout  (input) 

=  torpedo  speed  made  good  during  search  (input) 


integral  part  of  [f^  (r 


c(o) 


Jo  = 


integral  part  of  [f^  (r 


i(o) 


-  d  )  /  (s  At)] 

8 

-  dg)  /  (s  At)] 


s . 
J 


j  £  j' 


j  >  j- 


For  superscript  j  =1,  2,  3,  ... 

N(j)  =  N(°)  +  jx  +  (j  -  !)  C 


A-l 


Arthur  I )  Little.  Inc 


For  subscript  j  =  1,2,...,  n(2)-  n(o) 


0.  =  0 
J  o 


XTj  XTj-l  -  sj  At  sin  Gj»  yTj  -  y.f.  2  +  s  At  cos  0 


Tj-1  cos 


0(o)  =  o(1)  =  0 


Ntaai  1?  C0"PUt‘!  the  "C°rrected"  syro  course  angles  0<2),  0«> 
0  *  iarge  loops,  lor  „hich  rhe  one  lot  .«  ls  J 

For  subscript  j  =  N(m)  -  N(o)  -  C  +  1,  N(m>  _  N(o) 

0  = 
j 

XTj  =  ^j-l  +  sj  sin  ®- 
yTj  =  yTj-l  +  sj  At  cos  ®- 


For  i  = 


3  !>•••»  max  j  compute 

xj  '  x3+n(0)  '  ?j+N<0)  •  y/  *  yJ+N<°)  -  nJ+1)(o) 

1  x 
xTj  =  xtj  ~Vn(o)  >  yTj  =  yTj  -  n.+Nco) 

r3  ’  1(X01)2  +  <^2)1/2  •  -  t<4/  +  <4>V 

sin  B  =  x1  /  r.,  cos  B  =  y1  /  r 
3  3  J  j  yj  j 

Sin  *«  ‘  4)  '  rTj  •  “s  BTj  *  yjj  /  rXj 

Find  B  and  B  ,  -tt<B  <tt  ,  -tj<b  \  <  n 
J  Tj  J-  Tj 
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Arthur  D  Little,  Inc 


* 


* 


sin  0o(m)  =  sin  0  ^  cos  -f  cos  0  ^ 

8  L 


cos  0g^  =  COS  ©L^  cos  -  sin 


Find  0  ^  in  radians,  -  t,  <  0  ^  <  n. 
8  S  - 


Compute 


sin 

sin 


a o(m)  =  0  (m)  - 
g 

=  proportion  in  Proportional  Navigation  (input) 

^©max  =  ^-im:i-tat;ion  on  correction  (input) 

0(m)  =  0(m-l)  +  min  (|A0(a)  |  AQ  }  if  A@(®)  >  0 

0(m)  =  0(m-l)  _  m.n  (|Ae(»)|  A0  }  -f  &sCa)  <  Q 

4  1  1  max 

This  completes  the  loop  to  find  0^°\  Starting  with  ©^  *  ©^  =  0  , 
repeat  the  loop  to  find  0^2\  0^\...,  0^aax 

4.  Remainder  of  Computations 

The  remainder  of  the  computations  are  made  fey  the  procedures 
described  in  Sections  E,4,  E.5,  and  E.6  of  Appendix  E  of 
Reference  (a) . 
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